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In this work, we construct a phenomenological effective model for the heavy-light quark systems, 
which consist of (u, d, c, b) quarks, i.e. extended nonlocal chiral-quark model (ExNLChQM) to com- 
pute the heavy-meson weak-decay constants fo and /s. ExNLChQM is based on the heavy-quark 
effective field theory as well as the dilute instanton-vacuum configuration. In ExNLChQM, a certain 
portion of the heavy-meson mass is considered to be generated from the nontrivial QCD vacuum 
contribution, similar to the light quarks in usual instanton approaches. Hence, the effective heavy- 
and light-quark masses become momentum-dependent and play the role of a smooth UV regulator. 
Employing a generic external-field method applied to the effective action from ExNLChQM, we 
obtain fo = (169.28 ~ 234.57) MeV and fs = (165.41 ~ 229.21) MeV from the numerical results, 
depending on different model parameters. These values are in relatively good agreement with ex- 
perimental data and various theoretical estimations. We also discuss the heavy-quark effects on the 
QCD vacuum, and the decay constants fo' and /_b» in terms of the heavy-quark spin symmetry. 
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I. INTRODUCTION 

For a couple of decades, the heavy-quark effective field theory (HQEFT) has been known to be very powerful 
theoretical tool to investigate various heavy mesons and baryons, including the charm and bottom quarks 
HQEFT has been well tested in comparison with experiments for those heavy particles by the various collaborations, 
such as D0 [4], BaBar [5], Belle [6], BES [7], and so on, and accumulated successful results. HQEFT are based on two 
heavy-quark symmetries, i.e. spin and heavy-flavor ones. In general, they are not obvious in full QCD and become 
intact at mq — >■ oo, where mq stands for the heavy-quark current mass. As easily understood, the theory contains 
a small expansion parameter l/mg so that QCD can be expanded systematically in terms of it. Additionally, one 
usually defines the heavy-quark velocity v, satisfying = 1, which designates the velocity superselection rule [T]. 
The heavy-quark field can be redefined with v by integrating out the irrelevant (small) component of the spinor. As 
a consequence, at the leading l/mg expansion, the free heavy-quark Lagrangian manifests the spin and heavy- flavor 
symmetries. At the hadron level for instance, due to the spin symmetry, there appears mass degeneracy for spin-0 
and spin-1 mesons with the same parity. Interestingly, only a few percent deviation from the mass degeneracy has 
been observed in the real heavy- meson spectra [8 . The heavy-flavor symmetry becomes obvious as well, since the mg 
dependences are removed from the heavy-quark interactions at v. 

Many effective approaches have been made so far to study the heavy-light quark systems, such as D and B mesons for 
instance, based on HQEFT [51116). On top of HQEFT, the light-quark components inside the heavy meson were treated 
in various ways. Since the light-quarks are believed to be governed dominantly by spontaneous breakdown of chiral 
symmetry (SBCS), the Nambu-Jona-Lasinio model (NJL) type interactions were employed frequently, i.e. manifesting 
the chiral symmetry and its dynamical breakdown [9]. Similarly, other nonperturbative QCD approaches were also 
taken into account to reserch the heavy-light quark systems: Bethe-Salpeter (BS) approach [TOIIII], relativistic quark 
model (RQM) |T2], field-correlator method [13;, light-cone formalism [T3], QCD sum rule (QCDSR) T5], covariant 
chiral-perturbation theory (ChPT) [16j, and so on. As the numerical evaluation of the first principles of QCD, a 
considerable number of lattice QCD (LQCD) simulations were performed for the related subjects as well [17H22]- 
A mesonic effective interactions were developed for the heavy and light mesons with the flavor symmetry-breaking 
effects 

Among the physical observables to be addressed using these effective and LQCD approaches for the heavy-light 
systems, the heavy-meson weak-decay constants for D and B mesons have been investigated intensively as in the 
above references, since they are related deeply to the CP violation [21] as well as the Cabbibo-Kobayashi-Maskawa 
(CKM) matrix [5S]. Moreover, the various decay modes for those heavy mesons, such as the (non)leptonic decays. 
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are of importance in understanding the heavy-light systems profoundly. We note that CLEO collaboration have done 
various experiments for these physical quantities |26H28) , and the recent experimental data from CLEO were analyzed 
carefully in Refs. [IHIEOIj giving an estimation = (206.7 ±8.9) MeV. Using the B~ t^D data from the Belle [5] 
and BarBar [S] collaborations, the value of /s was estimated by (204 ± 31) MeV [13]. 

In the present work, considering the increasing interest on the heavy-light quark systems, we want to construct 
a new effective model for this subject. We note that the model is based on two ingredients in principle: HQEFT 
and nonlocal chiral-quark model (NLChQM) [5TI - I55] for the heavy and light quarks, respectively. Since HQEFT 
has been explained briefly already, we focus on NLChQM here. It is worth mentioning that NLChQM defined in 
Minkowski space is derived from the instanton-vacuum configuration. Instantons have been believed to be one of 
the important agents responsible for SBCS in the low-energy region [39l - l4T| . In the instanton model, the light quark 
acquires its effective mass ^ 350 MeV, which is a consequence of the nontrivial interactions between the quarks and 
(anti)instantons, as an indication of SBCS as well. Those nontrivial interactions also give the momentum-dependency 
to the effective-quark mass, which plays the role of a natural UV regulator. There are only two intrinsic parameters 
in the instanton model: average (anti)instanton size p « 1/3 fm and inter-(anti)instanton distance i? « 1 fm [39j . 
As long as we are staying at the SU(2) light-flavor sector for (u, d) quarks, these parameters are well determined 
phenomenologically and can reproduce various nonperturbative-QCD quantities qualitatively very well |31H33j . It 
is worth mentioning that one of the distinctive features in the instanton model is the interactions between the light 
quarks and pseudoscalar (PS) meson, which takes part in the model by an appropriate SU(2) bosonization process [42], 
becomes nonlocal ^43J. Hence, NLChQM derived from the instanton model bears this distinctive feature naturally. 

Keeping these ingredients in mind, we pay attention to the the mass difference between the heavy meson and heavy 
quark. In the limit of the infinite heavy-quark current mass, vtlq — )■ c», these two masses are the same. However, 
in reality, one observes they are differ from each other [8|. Intuitively, the heavy-meson mass can be regarded as the 
total sum of the various light and heavy quark masses, i.e. the effective and current ones. Here, we assumed the chiral 
limit so that the light-quark current mass is neglected. Below is an example for D-meson mass: 

1869MeV > 350 MeV -I- 1290 MeV . (1) 

D meson mass light quark effective mass lieavy quark current mass 

From Eq. ([!]), we see that about 200 MeV mass must be generated additionally, even there is a certain value for the 
binding energy for the heavy meson. In this sense, we assume that this additionally necessary [effective) mass comes 
from the similar mechanism with the light quark, i.e. from the nontrivial contributions from the QCD vacuum. If 
this is the case, we take into account that the effective heavy-quark mass can possess a similar role as a UV regulator 
as that for the light quark. Considering all these ingredients, we write an effective interactions (action) , based on the 
ideas of HQEFT and NLChQM, extended to the heavy-light quark systems, i.e. extended NLChQM (ExNLChQM). 

In order to test the validity of the present model, we compute the D- and _B-meson weak-decay constants, taking 
into account their essential implications in physics as mentioned above, using ExNLChQM. We choose all the model 
parameters from basic NLChQM for the light quarks which are the values for the renormalization scale and light-quark 
effective-quark mass at = 0, where q denotes the transferred momentum ^43j . As for the momentum dependences 
of this mass, we try several possible different choices, as will be discussed in Section II. An unknown parameter from 
the heavy-quark side, i.e. the effective heavy-quark mass at = is simply determined to satisfy Eq. |l| with small 
light-quark current mass. As noticed, since the effective masses are decreasing functions and play the roles of UV 
regulators, there is no need to include any adjustable free parameters by hand, such as the three-dimensional cutoffs. 

Once we obtain the effective action of ExNLChQM, it is an easy task to evaluate the decay constants via conventional 
functional treatments. From the numerical calculations, we obtain — (169.28 ^ 234.57) MeV and /s = (165.41 ^ 
229.21) MeV, depending on different model parameters for the momentum dependences of the effective heavy and 



light quark masses as shown in Eq. ( 10 ) in the next Section. These values are in relatively good agreement with 
experimental data and various theoretical estimations. By increasing the renormalization scale for the heavy quark, 
we can obtain more compatible numerical results with the experimental data, which indicates that the heavy quarks 
make effects on the QCD vacuum to a certain extent, different from the light quark. We also discuss the vector 
heavy-meson weak-decay constants fD* and /s* in terms of the heavy-quark spin symmetry. It turns out that the 
ratios, fD*,B-'/fD,B are in 1.01 ~ 1.03, which are relatively smaller than those from other theoretical calculations. 

The present work is structured as follows: In Section II, we discuss how to formulate the effective model for the 
present purpose with brief explanations on HQEFT and NLChQM. Here, we also elucidate the parameters for the 
present model. Section III is devoted to the numerical results with relevant discussions. We close the present work 
with summary and future perspectives in the final Section. 
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II. GENERAL FORMALISM 



In this Section, we want to explain the motivation for ExNLChQM and how to construct it based on HQEFT and 
instanton physics. First, we briefly introduce HQEFT following Refs. [1]. The heavy-quark symmetry is fully satisfied 
in the limit of mg — > oo, where mg stands for the current-quark mass for the heavy-flavor SU(2) quarks Q = {c,b). 
Since the top quark, whose mass is about a hundred times larger than these two quarks, we exclude it from the present 
discussions. Conventionally, the heavy quark filed Q can be rewritten as a soft one Qy with a definite velocity v [IK3]: 



imnv-x 



Qv (^) ; 



(2) 



where we define the light-like vector v satisfying = 1. Moreover, we write v — (1, 0, 0, 0) for definiteness, indicating 
that the heavy quark is at rest. Using this definition of the heavy-quark field, the Dirac equation for it can be modified 
as [J 



Q (ip ~ mq) Q = Qv{iv ■ d)Qy. 



(3) 



Taking into account the spin symmetry of HQEFT, the pseudoscalar (scalar) and vector (axial- vector) mesons are 
grouped as a spin doublet H {K) [5]: 



H 



K : 



1+^ 



(4) 



Here, {S,V,V'^,A'^) correspond to the fields of the (scalar, pseudoscalar, vector, axial-vector) mesons, consisting of 
the quark contents qQ minimally. In the present work, we will work with the SU(2) light-flavor sector, i.e. q = (u,d). 
Since we are interested in the weak-decay constants for the PS mesons, such as D and B, in order to test the present 
model, we will concentrate on the doublet field H only hereafter. 

Now, we take a look on the light-flavor sector. It has been well known that the light flavors are strongly gov- 
erned by the spontaneous breakdown of chiral symmetry (SBCS). Among the various chiral effective models, the 
instanton-inspired models have accumulated successful achievements to investigate important nonperturbative QCD 
properties and hadrons (39j Je!1. ,43J. Based on the dilute-instanton model, the nonlocal chiral quark model (NLChQM) 
was developed and has been applied to many nonperturbative problems [3Tti37l [37l l38] . Considering its successful 
applications so far, we will employ NLChQM for the light-flavor sectors. The effective chiral Lagrangian of the model 
has the form as follows: 



off 



(5) 



where rhq = diag(m„, m^) and Mq denote the current-quark and effective masses for the light quarks. According 
to the isospin symmetry, we choose — md = 5 MeV throughout the present work [S]. Note we do not consider 
the strange quark here whose mass is about (100 ~ 200) MeV [8], so that the weak-decay constants for Ds and 
Bg will not be taken into account. The reason for the exclusion of the strange quark is as follows: The inclusion 
of the strange quark in Eq. ([5| breaks the validity of the effective Lagrangian in the leading Nc contribution |44j . 
To remedy this problem, one needs to go over the leading Nc contribution, and it turns out that the meson-loop 
corrections, corresponding to the large- iVc corrections, are necessary to remedy the problem |44fl46j . Hence, we will 
confine ourselves to the SU(2) light-flavor sector in the present work to avoid theoretical complexities. Note that 
Mq is generated from the nontrivial quark-(anti)instanton interactions in terms of the instanton physics, manifesting 
SBCS as well as the QCD- vacuum contributions 43J: 



Mq - M^F^it), Fit) = 2t 



Ut)Ki{t) - hit)Koit) - -h{t)Kiit) 



(6) 



We have used the notation t — \p\p/2, and la and Ka stand for the modified Bessel functions with the order a. Note 
that C/5 plays the role of nonlinear PS-meson fields and reads: 

2F3 



C/5 = exp 



»75(t 
F^ 



= 1 



ilbj-r ■ tt) 
F^ 



(7) 



Here, F^^ indicates the light PS-meson weak-decay constant, whose empirical value amounts to F^^ « 92.3 MeV. As 
for the interaction Lagrangian in Eq. ([5|, we want to emphasize on that the quark and PS meson interactions are 
momentum-dependent, i.e. nonlocal interaction, represented by the third term in the square bracket in the right- 
hand-side of Eq. ([5]). 

Taking into account all the ingredients discussed for the light and heavy sectors so far, now, we are in a position to 
extend this NLChQM to a system with the heavy-light quarks. Our strategy for this purpose is as follows: 
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• We extend the flavor SU(2) symmetries for {u,d) and {c,b) into that for SU(4), consisting of {u,d,c,b) quarks. 

• The form of the interaction in Eq. ([s]) is assumed to be kept even for the flavor SU(4), manifesting momentum- 
dependent quark-PS meson couphngs. 

• Hence, the corresponding heavy PS mesons are represented by the nonhnear form as in Eq. 0. 

• We consider that the heavy quarks also possess their effective mass as the hght one does. We wifl discuss this 
in detail soon. 



As has been known well, the flavor SU(4) symmetry is strongly broken in reality, mu,d ^ "^c,b- For instance, in 
Ref. [23], the mesonic effective Lagrangian was constructed in terms of the SU(4) symmetry with the {u,d,s,c) 
flavors, and the symmetry-breaking effects were taken into account by suppressing the virtual heavy-meson exchanges 
DC m0/m$, where and m$ denote the light- and heavy-meson masses, respectively. Being different from Ref. 
however, we are treating the quark degrees of freedom (d.o.f.) in the present work, and the symmetry breaking 
effects are included explicitly and microscopically by the current and effective masses for the heavy and light quarks, 
although the current heavy-quark mass will be subtracted from the resultant effective Lagrangian by the heavy-quark 
symmetry [T]. Note that the physical heavy- meson mass is used to compute the effective heavy-quark mass as in 
Eq. (19) as well as the numerical input as in Eq. (23), manifesting the symmetry-breaking effects. Moreover, the 



symmetry breaking effects in the present framework can be deduced from the different UV cutoff masses for the heavy 
and light sectors, as will be discussed in Section III. 

Upon this strategy, we can write an effective Lagrangian density for ExNLCHQM in a neat form as follows: 



^ExNLChQM _ ^ 



(8) 



where ip is for the SU(4) flavor spinor, i.e. = {u, d, c, 



while Dp for the covariant derivative iD^ = 



id,,. 



Here, we use the notations Cq^q for electric charge for a quark {q, Q) and for U(l) gauge held. Taking into account 
the SU(2) isospin and SU(2) flavor symmetries for the light- and heavy-quark sectors [T], we can simply rewrite the 
spinor for a two-component column matrix as 



tfj = {u,d,b,cf {q,Qf 



(9) 



The current-quark mass matrices read "rhq — diag(TO„, m^, 0, 0) and riiQ = diag(0, 0, rric, mi,) again. The square-root of 
the effective quark-mass matrix \/A^ is also defined by VA4 — diag{y/M^, y/Md, y/Mc, y/Mh), in which Mf^^ Q-^ stand 
for the effective (light, heavy)-quark mass. In the present work, instead of using Eq. (|6]), we employ parameterized 
forms for the effective masses for analytical and numerical convenience as follows: 



(9,Q) 



(9,Q),0 



A2 



Hq-.Q) 



A2 



^(<!,Q) 



(10) 



where Mq corresponds to the soft-component of the heavy-quark field Qy as understood. Note that, from now on, we 
choose TLq^Q as positive-integer free parameters. As a trial, we choose the cases for (n^, uq) = 1 and (n^, uq) — 2, being 
assigned as Model I and Model II, respectively. As mentioned in the previous Section, A denotes the renormalization 
scale of the present model and relates to the inverse of the average (anti)instanton size p ~ 1/3 fm, resulting in 
A « 591 MeV. In Figure fl] we depict the effective light-quark mass in Eq. ( [lO| ) as a function of the momentum 
transfer q for n — 1 (solid j7 2 (dot), and 3 (dash) for Mq g = 345 MeV in Euclidean space [33]. We also draw that 
from the dilute-instanton model (long dash), given in Eq. ^6]). As shown in the figure, when n = 2, the parameterized 
mass is very similar to that from the instanton model. Moreover, as n increases, the mass curves get decreasing faster. 
As Model III, we set nq = 2 for the light quark and hq = 1 for the heavy quark, i.e. {nq,nQ) — (2, 1), considering a 
possibility that the momentum dependences for the heavy and light quarks may different from each other. We will 
not take into account the case with {nq, uq) = (1, 2), since the numerical results are almost the same with those from 
Model HI. The reason for this similarity can be easily understood by seeing Eqs. (24) and (25) in Section HI, since 
the term (MqMq)^^'^ in the numerator governs the difference between Models. 
The nonlinear PS-meson field Uc, in a (4 x 4) matrix form reads: 



exp 



exp 

't75T'fl 

-Fn 



z75T-7r 



exp 
l2x2 exp 



t75T-n 

-Fn 



n' 



L2X2 + — — 

— 12X2 



(11) 



where -F,r,n,n' indicate the weak-decay constants for the isovector PS mesons (tt, H, H'), consisting of {qq, qQ or Qq, QQ) 
quark contents. Note that we have off-diagonal terms in which are necessary for the heavy-light quark mixing 
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FIG. 1: Effective light-quark mass Mq as a function of the momentum transfer q in Eq. (10 1 for n — 1 (solid), 2 
and 3 (dash) in Euclidean space. We also draw that from the dilute instanton model (long dash) in Eq. Q. Here, i 
A=l/p^ 591 MeV and M,,o « 345 MeV. 



(dot). 



terms. The subtraction in the off-diagonal terms, i.e. — l2x2, indicates that there is no flavor mixing between heavy 
and light quarks. We expand up to ^^(Tr, 11, 11'), since we will compute the matrix element (0|Jvk| PS meson) for 
Fd,b- Moreover, because we are not in interested in the quarkonia states such as J/'^ ~ cc and T ^ bh, which are 
about twice-times heavier than D and B mesons, we set the lower-right (2 x 2) block matrix to be a unit matrix, 
resulting in a Lagrangian for the heavy quarks in Eq. ([2]). However, it is a straightforward task to include those 
quarkonia states in the present framework. The explicit expressions for the (2 x 2) block flavor matrices in lA^ are 
written with the light and heavy PS mesons as follows: 



^J2^:- -tt" 



n = \/2 



5° B+ 

D- 



n = \/2 



D+ 

B- 5° 



(12) 



Also, the explicit expression for lA^ is given in Appendix in a (4 x 4) matrix form. Note that, if we replace the present 
(m, d, c, 6) flavors into (w, d, s, c), one recovers the SU(4) PS-meson fields given in Ref. [23]. It is worth mentioning that 
Eq. ([8| breaks chiral symmetry explicitly, whereas it is dynamically broken by the finite Mq^Q values, which indicate 
the nontrivial QCD vacuum effects [33 SOI HI] ■ This chiral-symmetry breaking pattern is inherited from NLChQM 
beyond the chiral limit [33] . 

Using the definitions from Eq. ^ to Eq. (12), the effective Lagrangian in Eq. ([8| can be represented in three 
separate parts, i.e. light-light (LL), heavy-heavy (HH), and heavy-light (HL, LH) quark terms: 



^ExNLChQM _ 



L 



LL 
cff 



HH 
cff 



c 



(HL,LH) 



q{ipq 
' 1 



cff 

Mq) q 



1 



m^/M~qq 



1\/M^ [75 (iT • tt)] 



HL 




Mq) Q] 



HH 



(13) 



LH 



where all the masses in Eq. (13) are 2x2 diagonal matrices for each flavor SU(2). Being similar to the heavy-quark 
field in Eq. ([2]), the zero-spin heavy PS-meson field can be redefined with its soft component as follows: 

n^, n = e""«^-^n„. (i4) 



n 



-imnv-x 



It is an easy task to redefine C^g in Eq. (13) using the Dirac equation for the heavy quarks in Eq. (|3|), resiflting in 



/'HH 

■'-efl - 



^{iv ■ d) - Mq] g, = [{tv ■ d) 

In the second step of Eq. (15), we use the projection relation ^Qy = 1^ 
^(HL,LH) ^j^^ following form: 



-Mq]Q,. (15) 
Using Eqs. ^ and (14 1, we can rewrite 



/'HL _ 
•'-off — ^" 



'Mqq, C]i^- 



Fu 



'M„ 



- l + i 

{iT ■ n„)— — 75 



(16) 
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which satisfies £^J^^'^^'> = C^J^'^^\ According to the spin symmetry for the heavy quarks and comparing with 
Eq. Q, the heavy PS meson can be re-defined in a isotriplet spin-doublet meson for (0~, 1~), as follows: 



1 + ^ 



75 (ir • n^,) 



1 + ^ 



T H. 



(17) 



Using a generic functional integral technique for the Grassmann variables given in Appendix, finally, we can arrive at 
an effective action for the heavy- light quark systems from the effective Lagrangian density in Eq. (13): 



S. 



LL+HL+LH 



-iSp In 



1 



1 



MqH^Mq {ivd-Mq)- 



H 



, (18) 



where Mq = niq + Mq, and Fh denotes the weak-decay constant for the spin-doublet mesons as in Eq. ( 17|. We have 
used notations tt = (t • tt) and H = t ■ H for convenience. It is worth noting that the effective action in Eq. ( 18 ) is 
in principle equivalent to the first term of Eq. (36) in Ref. [9 , except for the momentum dependent quark-PS meson 
coupling strengths. 



As a next step, we estimate Mq^q in Eq. ( 10 1 from a simple phcnomcnological analysis. In this consideration, the 
heavy PS-meson mass can be understood as 



M 



H 



[niq + Mq,o]l^ + [toq + Mq,o]h : 



(19) 



where we have ignored the binding energy for the meson. From the experimental data for D and B mesons, we can 
write 



Md = 1869.57 MeV « (7nc 
Mb = 5279. 17 MeV w {nib 



Mq,o + Mqa) = 1625.0 MeV + Mq^q Mq^q 
Mqfi + Mqa) = 5025.0 MeV + Mq^q ~^ Mq^ 



244.57 MeV 
254. 17 MeV. 



(20) 



The numerical inputs, which are taken from Ref. [5], are summarized in Table [ill shown in Eq. (20), it is necessary 
to add the effective heavy-quark mass to reproduce the heavy-meson mass appropriately. To estimate the values of 
Mqa above, we choose Mq^ « 350 MeV as a trial. If we think the binding energy for the mesons, the estimated value 
Mq q = (240 ~ 250) MeV must be the lower bound of its real one. From these observations and previous discussions, 
we assume the following two: 

• The finite effective heavy-quark mass, Mq q is generated from the same mechanism with that for the light 
quark-instanton interaction, representing the nontrivial QCD vacuum effects. For the numerical calculations, 
Mq^o is estimated by Eq. (19 1. 



The instanton ensemble is not affected much by the heavy sources such as the heavy quarks, resulting in that 
one can use Aw l/p~ 600 MeV again for Eq. (10). However, we will also discuss the possibility for A for the 



heavy quark to be changed from about 600 MeV in Section III. 



Hence, we use the effective heavy-quark mass in the momentum space as in Eq. (10), although we do not go over its 
microscopic derivation, which must be beyond the scope of the present work. 



rric 


Md 




Mb 


(1290.01E; °5) MeV 


(1869.57 ±0.16) MeV 


(4670.0tE!;J^) MeV 


(5279.17 ±0.29) MeV 



TABLE I: Numerical input values for the heavy quarks and mesons for Fd.b in Eq. (24 1, taken from Ref. 
the masses for the charged heavy mesons. 



Here, we choose 



III. NUMERICAL RESULTS AND DISCUSSIONS 



In this Section, we present numerical results, computed using the effective action of ExNLChQM in Eq. (18), to 
verify the validity of the present model. For this purpose, we perform simple calculations for the heavy PS-meson 
decay constants Fd^b through a conventional functional method. First, we define the (_B, D) = H meson weak-decay 
constants as follows |9|: 

(0|g-(a;)7^(l - i,)t''QAx)\H\p)) = 2ip^Fh5''\ (21) 



7 



Q(k) 




q(k-p) 

FIG. 2: Feynman diagram for the weak decay for a heavy meson H — {D, B) with the momentum p. Thick and thin sohd lines 
indicate the heavy Q = (c, b) and hght q = [u, d) quarks, respectively. F denotes an arbitrary Lorentz vertex depending on the 
decay processes. 



where the superscripts (a, b) indicate the isospin indices, while p is the on-mass shell momentum of H with the velocity 
V. To evaluate the matrix elements in Eq. (21 ), it is convenient to employ the external-current method, which can be 



represented by the following effective action: 



CLL-I-HL+LH 



ip - Mq 




MqH^Mq + an ■ Jh ] (iv ■ d - M, 



(22) 



Here, a and J" denote the relevant external source and current with a certain isospin index a, respectively, in general. 
Although we do not evaluate F^^ numerically using the effective action in Eq. (22 1 in the present work, since we have 
done it several times in the previous works [47] , the numerical result for turned out to be about 93 MeV, which 
is in a very good agreement with its empirical value, with the phenomenological instanton parameters A = 591 MeV 
and Mq^o = 345 MeV. These values wiU be employed further for Fd,b- 

To evaluate Fd,b, one needs to determine the value for v ■ p, considering the definitions in Eq. (21 1. From the 
on-mass shell condition for the meson H at v [H], one is lead to 



■ p « AMd,b = Md,b - rric^b- 



(23) 



We note that this quantity in Eq. ( 23 ) describes the long-distance nonperturbative effects 
quark rest frame v^i = (1, 0, 0, 0), one can write Pf_i = (AMg, 0, 0, 0) 



. Hence, taking the heavy- 
Performing standard functional derivatives with 



respect to an and H for the weak current in Eq. (21 ), the analytic expression for Fh can be obtained as follows: 



F' 



H 



4A^ 
V ■ p 



d'^k yjMq^k-pMq^k v ■ k-v ■ p 
(27r)4 [(fc-p)2-M2,_p]^ 



M, 



(24) 



Note that, in the right-hand-side of Eq. (24), we replace the effective heavy-quark mass into that at fc^ = in the 
denominator, i.e. Mq — >■ Mq^, since this replacement does not make any considerable changes, while the numerical 
calculations get easier to a great extent. It is worth mentioning that, due to the decreasing behavior of the effective 
quark masses in Eq. (10 1 with respect to momentum transfer, the integration over k in Eq. (24 1 does not diverge 



even without any cutoffs included by hand. The integration in Eq. (24) can be evaluated by 1) performing the Wick 
rotation to Euclidean space and 2) Cauchy integral over k'^ = k^ — Mq,o, as done in Ref. [9j. After these analytic 
calculations, we arrive at an expression for the heavy-meson weak-decay constant: 



AM, 



H 



2]2n„ 



k^dk y/Mq^oMQ^o{AMH - MQ,o)[fc^ + (Mq,o - AMg)^ + n^A^] 

^2 [fc2 + (Mq,o - AAfe)2 + n,A2]2»,[fc2 + ^Mq^o - AMh)2] + (71,A2)2«QM2^ 



(n,A2)"«/2 



[fc2 + (AfQ,o - AAfo)2 



jA2]"./2 



(nQA2)"o/2 

(fe2 nQA2)»Q/2 ' 



(25) 



where «<j,q are defined with Eq. (10). 

In order to compare the present numerical results with other theoretical ones, we normalize them alternatively as 
Fh Jh = V'2.Fh- Then we have the following numbers for the heavy-meson weak decay constants, fD,B, for Model 
I with {nq^nq) = 1: 



Jd = 234.57 MeV, /s = 229.21 MeV. 



(26) 
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FIG. 3: fo (left) and fs (right) computed by various theoretical approaches, which are indicated on x axis and given in 
Table The shaded areas stand for the experimental data with corresponding errors [141 1291 [30] . The horizontal thin-solid 
lines indicate the center value of the experimental data. 




FIG. 4: fo (left) and fs (right) for different renormalization scale A in Eq. ( |10[ ) for the heavy quark. Here, we fix it for the 
light quark as A = 591 ~ 1/p MeV. The shaded areas stand for the experimental data with corresponding errors [141 1291 [30] . 
The horizontal thin-solid lines indicate the center value of the experimental data. 



Similarly, those values computed via Model II with (uq, uq) ~ 2 read as follows: 

/d = 169.28 MeV, /s = 165.41 MeV, (27) 
and, from Model III with [nq^nq) ~ (2, 1), we have 

/d = 189.53 MeV, /s = 185.20 MeV. (28) 

Note that all the numerical results of Model I, II, and III are evaluated at A = 591 MeV for the both heavy and 
light quarks. If we take into account the errors in heavy-quark and heavy-meson masses |8j, these values would be 
changed by ±(2 ^ 3) MeV. In Table [TTj we list experimental and other theoretical results in comparison with ours. 
For this table, we refer mainly to Ref. |14) . From the table we see that our Model-I results are in about 10% deviations 
from the experimental data f(D,B) = (206 ± 8.9, 204 ± 31) MeV [29l [3^, whereas those from Model II considerably 
underestimates the experimental data by about 20%. Model III results locate in the middle of those from Model I and 
II, as expected, and closer to the experimental data than others. The theoretical results from the lattice QCD (LQCD) 
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with the improved clover action in the quenched approximation, they had f(D,B) = (211 ± 14, 179.0 ± 18) MeV |19) . 
A field-correlator method (FC) was also employed then gave f[D,B) = (210 ± 10, 182 ± 8) MeV [H]. Comparing with 
these values, our results are relatively compatible with them, showing only several percent differences. The QCD sum 
rule (QCDSR), J^d^b) = (195±2G, 206± 15) MeV P^, and Bethe-Salpeter method (BS), f(D,B) = (230±25, 196±29) 
MeV [inilll], also present the comparable results with ours. Note that there is tendency Jd < /b from the QCDSR 
results, being different from others. The relativistic quark model (RQM) [12] and extended NJL model (ExNJL) [9] 
present relatively large values for fjj in comparison to other estimations. Note that the theoretical framework of 
Rcf. 9J is in principle the same with ours, except for the quark-meson coupling schemes. Although we did not list in 
Table there are many other LQCD results, such as fs = (204 ± 8) MeV 0] and (216 ± 9) MeV [22, which locate 
inside the experimental errors. 

In order to test the A dependence for the heavy quark, considering the possibility for the QCD vacuum to be 
affected by the heavy sources differently from the light quarks, in Figure |4j we plot the curves for fu (left) and fs 
(right) as functions of the renormalization scale A in Eq. ( 10 ) for the heavy quark. Here, the renormalization scale for 
the light quark remains unchanged, i.e. A « 591 MeV, since it has been well known from the usual instanton models 
theoretically as well as phenomenologically [321 SHI 133 . In general, as A increases, the values for fo.B also does 
smoothly. As shown in the figure, Model I results are almost excluded from the experimental data for any of A values. 
On the contrary, if A becomes about 1.0 (0.8) GeV for Model II (III), one obtains considerably good agreement with 
the data. From this observation, one can argue that the renormalization scale for the heavy quark can be modified by 
the interaction between the heavy quarks and QCD vacuum into a larger one, i.e. A ~ 0.6 GeV — 1 GeV, which may 
indicate the effects of the heavy sources to the vacuum. Since A is proportional to the inverse size of (anti)instanton, 
the increasing of A stands for the spatial shrinkage of the (anti)instantons. We notice that this tendency is consistent 
with that the effective heavy-quark mass ^ 250 MeV, estimated by Eq. pO| ), is smaller than the effective light- 
quark mass ~ 350 MeV, since the absolute value for the effective mass is proportional to the interaction strength 
between the quarks and (anti)instantons. As understood by seeing Figure |4j if we tune the A values for the heavy 
quark, we can reproduce fo.B being compatible very well with the center values of the data. For instance, we obtain 
f{D,B) = (205.00, 200.32) MeV with A = 1.0 GeV for the heavy quark via Model II. Note that these values show 
only a few percent deviations from the experimental data. However, we do not provide those numerical results in the 
present work, because we are interested in showing the validity of ExNLChQM, not in those fine tunings. We want 
to leave such a phenomenological studies for the heavy- light systems for the future works. 

By construction of the present model with the spin symmetry, the decay constants for the vector heavy mesons 
/d. and fs* can be computed using Eq. ( 25 ) as well. Only difference between the pseudoscalar and vector mesons in 
the present theoretical framework is their masses as the inputs. Here are some discussions for the vector mesons: We 
note that the vector-meson d.o.f. are not incorporated with generic NLChQM, since all other fields except for the PS 
one are integrated out [H] , resulting in that the present ExNLChQ VI is not accompanied with them by construction. 
There have been several attempts to take into account the vector mesons within NLChQM in terms of the vector- 
meson dominance (VMD) [48]. In that way, ExNLChQM can be modified for D* and B* without considering the 
spin doublets in HQEFT. One can also include the heavy vector- meson fields nonlinearly in the present model as 
in Eq. (11), being similar to Ref. [23]. In order to test the VMD-modified ExNLChQM, one may compute fD*,B' 
as well as the heavy vector- meson light-cone wave functions [T3]. We would like to leave this VMD modification of 
ExNLChQM as future works. 

As the numerical results, ignoring the errors and uncertainties in the quark mass values, we have fw — 
(240.72, 172.74, 194.52) MeV and fo* = (231.80, 167.31, 187.31) MeV for Model (I, II, III). These computed val- 
ues for the vector mesons are also listed in Table [TTj with comparison with other theoretical estimations. In Ref. [14j . 
the values for f(D,B)* were computed by using the experimental values for f(D,B), which are underlined in the table, 
as inputs. Note that, from our model, the values for f(D,BY are about a few percent larger than those for the PS 
mesons. Similar tendency, f(D,B)* > f(D,B)i is also observed in other theoretical calculations as in the table. The 
ratio for the PS and vector heavy-meson decay constants, f{D.B)*lf{D,B) is a good quantity to see this tendency 
quantitatively. As seen in the table, other theoretical calculations provide the ratio, which ranges from 1.1 to 1.48 
with about 10% uncertainties. Being different from them, our numerical results show small numbers for the ratios as 
explained above. This smallness may indicate again the necessity to modify A for the heavy quark. 

From these observations discussed above, we can conclude that the present model, i.e. ExNLChQM reproduces the 
weak-decay constants qualitatively, in comparison to various experimental and theoretical results. Moreover, we 
have shown that there is a more room to improve the model from a phenomenological point of view by considering 
the heavy-quark effects on the nontrivial vacuum. 
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[MeV] 


Model I 


Model II 


Model III 


LCF [H] 


LQCD HH 


FC [13] 


QCDSR [IS] 


BS [lOldT] 


RQM [12] 


ExNJL [9] 


Id 


234.57 


169.28 


189.53 


206 ±8.9 fM 


211 ± 14 


210 ± 10 


195 ± 20 


230 ± 25 


234 


300 


fB 


229.21 


165.41 


185.20 


204 ±31 (30, 


179 ± 18 


182 ±8 


206 ± 20 


196 ± 29 


189 


180 


fo- 


240.72 


172.74 


194.52 


259.6 ± 14.6 


245 ± 20 


273 ± 13 




340 ± 23 


268 




fB* 


231.80 


167.31 


187.31 


225 ± 38 


196 ± 24 


200 ± 10 




238 ± 18 


219 




fo'/fo 


1.03 


1.02 


1.03 


1.26 ± 0.12 


1.16 ±0.14 


1.30 ±0.13 




1.48 ±0.09 


1.15 




fB'/fB 


1.01 


1.01 


1.01 


1.10± 0.16 


1.09 ±0.13 


1.10 ±0.13 




1.21 ±0.06 


1.16 





TABLE II: The values for fo.B and fD*,B* from experimental [291 130] (underlined and inputs for Ref. [14]) and theoretical 
results [MeV]. We also summarize the ratios fD*,B'/fD,B for various cases. All the numerical results of Model I, II, and III 
are evaluated at A = 0.591 MeV for the both heavy and light quarks. The explanations for the abbreviations and numerics are 
given in the text. 



IV. SUMMARY AND OUTLOOK 

In this vifork, we have constructed a phenomenological heavy-light quark effective model, i.e. ExNLChQM and 
tested it to reproduce the heavy-meson weak-decay constants. The present model was based on and motivated by 
the heavy-quark effective field theory and instanton-vacuum configuration. In what follows, we list some important 
points of ExNLChQM and observations of the present work: 

• We focused on the SU(4) heavy-light flavor sector, i.e. {u, d, c, b) quarks, since the inclusion of the strange 
quark breaks down the validity of the usual leading- iVc instanton-induced effective model and the top quark is 
too heavy to be incorporated with these four quarks. The flavor SU(4) symmetry-breaking effects are taken 
into account explicitly by the current and effective masses for the heavy and light quarks, using the physical 
heavy-meson masses as numerical inputs. 

• We assumed that the heavy quark also possesses its effective mass, generated from the interactions with the 
nontrivial QCD vacuum, being similar to the light quarks, resulting in the momentum-dependent effective heavy- 
quark mass, which plays the role of a UV regulator. As a consequence in the present model, the heavy meson 
and heavy-light quarks are coupled nonlocally. 

• The renormalization scale of the model was set to be A « 600 MeV, which relates to the inverse of the average 
(anti)instanton size. It was also approximated that the A for the heavy-quark effective mass is the same with 
that for the light quark, taking into account that the QCD-vacuum structure is not affected much by the heavy 
quarks, although the heavy sources may distort the vacuum to a certain extent. 

• On top of the heavy-quark spin symmetry for the heavy mesons, the D- and i3-meson weak-decay constants are 
computed using the ExNLChQM effective action, resulting in fo = (169.28 234.57) MeV and fs = (165.41 ~ 
229.21) MeV from the numerical results, depending on the instanton parameters. These values are compatible 
with experimental and other theoretical estimations. 

• As the renormalization scale A for the heavy quark gets increased, one obtains more compatible results with 
the data, although we do not provide those numerical numbers. This tendency indicates that the QCD vacuum 
can be affected by the heavy sources with different strengths in comparison to the light-quark case. 

• The vector heavy-meson weak-decay constants fu* and /s* are also computed in terms of the heavy-quark spin 
symmetry. From our calculations, the ratios fo* / fo and fs" / fs become only a few 1.01 ~ 1.03, which is rather 
smaller than other theoretical calculations, manifesting the spin symmetry. 

As a result, we can conclude that the present phenomenological model, ExNLChQM is quite reliable and promising to 
investigate the heavy-light quark systems, although there are several theoretical assumptions and improvements to be 
addressed from a microscopic point of view. For more detailed tests for the present model, theoretical computations 
for the decay constants for strange heavy mesons and various structure functions, such as the Isgur-Wise function, 
are under progress, and the results will appear elsewhere. 
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appendix 



The explicit expression for in the right-hand-side of Eq. (11) is given in (4 x 4) matrix as foUows: 



Us 



1 + 

^ qq 


i\/275 7r + 




qq 


FqQ 




F - 

r qq 




r qq 


FqQ 






1 


FQq 


FQq 




i\/2'i5B" 





FQq 


FQq 



FqQ 
FqQ 



1 



The functional integral with two Grassmann variables can be performed as follows: 

DQ,DQ,DqDqe'f'"-^-^'i+'^^''Q"+'^^^'i+'^^'Q^ = J DQ,DQ,DqDqe' f ''"-^'^+^^'''''^^'^'^+'''^^^+^^^-^''''^^'' 
= Det[6] / DqDqe'^'^'^'i^^-'^''^'^^i = Sp ln[S] + Sp ln[yl - C^S+c]. 
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